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We consider two systems A and B that share Einstein-Podolsky-Rosen (EPR) steering correlations 
and study how these correlations will decay, when each of the systems are independently coupled 
to a reservoir. EPR. steering is a directional form of entanglement, and the measure of steering can 
change depending on whether the system A is steered by B , or vice versa. First, we examine the 
decay of the steering correlations of the two-mode squeezed state. We find that if the system B is 
coupled to a reservoir, then the decoherence of the steering of A by B is particularly marked, to the 
extent that there is a sudden death of steering after a finite time. We find a different directional 
effect, if the reservoirs are thermally excited. Second, we study the decoherence of the steering of a 
Schrodinger cat state, modelled as the entangled state of a spin and harmonic oscillator, when the 
macroscopic system (the cat) is coupled to a reservoir. 
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I. INTRODUCTION 

The sort of nonlocality we call “Einstein-Podolsky- 
Rosen-steering” mm originated in 1935 with the 
Einstein-Podolsky-Rosen (EPR) paradox f?j. The EPR 
paradox is the argument that was put forward by 
Einstein, Podolsky and Rosen for the incompleteness 
of quantum mechanics. The argument was based on 
premises (sometimes called Local Realism or in Einstein’s 
language, no “spooky action-at-a-distance”) that were not 
restricted to classical mechanics, but were thought essen¬ 
tial to any physical theory (5] . The EPR argument re¬ 
veals the inconsistency between Local Realism and the 
completeness of quantum mechanics. It does not in itself 
rule out all completions of quantum mechanics, that are 
compatible with local realism. Nowadays, after the work 
of Bell, we realise this ruling out of all local realistic the¬ 
ories can be done in some special experimental situations 
EM. Any realisation of the EPR paradox HM3 , as a 
special case of EPR steering, is nonetheless important in 
giving a concrete intermediate result: the fact that quan¬ 
tum mechanics without completion cannot be consistent 
with local realism. A great advantage of EPR-steering 
over Bell-type tests is that they are more accessible to 
mesoscopic or macroscopic systems. EPR steering tests 
have been quantitatively studied or proposed for opti¬ 
cal down conversion EMU, optical systems in nonlinear 
regimes near or at critical points es, atomic gas ensem¬ 
bles at room temperature 1 191 f2T| , Bose-Einstein conden¬ 
sates (BEC) fTTHST] . and opto-mechanics |32H35| . This 
may lead to experimental tests of Schrodinger cat-type 
states. 

The EPR paradox manifests as a strong correlation 
between the positions and momenta of two spatially sep¬ 
arated particles (or some equivalent correlations). These 
correlations, once one assumes Local Realism, become 
inconsistent with the uncertainty principle. Given the 
strangeness of the EPR correlations, a likely hypothe¬ 
sis would be that they do not exist in real physical sys¬ 
tems. Indeed, the sort of correlations needed for EPR are 
not easily realisable in experiment. To date, only a few 



Figure 1. Our question is this: We consider two well-separated 
systems, A and B, that possess the type of quantum corre¬ 
lation we call EPR-steering. Each system is independently 
coupled to a thermal reservoir, at a time t = 0. This coupling 
induces a decay of the EPR-steering. We ask how does the 
decay of the steering depend on: the coupling j a , 7 {, to each 
reservoir; the thermal excitation n a , rib of each reservoir; and 
the measure that gives the strength of the EPR-steering. We 
will consider two types of EPR. systems: The first is a two¬ 
mode squeezed state. The second is a mesoscopic/ macro¬ 
scopic entangled “Schrodinger cat” state. 


experiments can justify the claim of loophole-free EPR 
correlations |36| . or steering without detection efficiency 
loopholes EH E3 [HD [35]. This paper examines why this 
is so: There are two possibilities. 

(1) It could be that quantum mechanics needs modi¬ 
fication to the extent that EPR correlations cannot be 
predicted. If we accept irrefutable evidence now exists 
for EPR correlations, then this probability would appear 
falsified. However, we comment that the experimental re¬ 
alisations for EPR steering have been for optical systems, 
not yet for atoms or mesoscopic devices. 

(2) Or else: a popular belief is that quantum mechan¬ 
ics is correct as it is, and predicts exactly why, to a great 
accuracy, the EPR correlations are extraordinarily diffi¬ 
cult to measure. In that case, these predictions need to 
be evaluated and verified by experiment. 
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This last point is what we examine in this paper: We 
study in detail why, according to quantum predictions, 
the EPR correlations are not easily realisable. As with 
the well known “Schrodinger cat” |.‘19lPIB] , we expect this 
is so, because of the effect that occurs when a quantum 
system is coupled to its environment, modelled as a large 
system - a reservoir [49] . This effect is called “decoher¬ 
ence”. 

Here, we examine some different sorts of decoherence 
that occur for EPR steering. As with the decay of en¬ 
tanglement, there are many cases one can study. It is 
worth noting that EPR steering was realised using high 
efficiency detection for optical amplitudes [14] [15], and 
has also been unambiguously detected in photonic sce¬ 
narios with exceptional losses (~ 87%) [SB]. Tests have 
also begun for its genuine multipartite form, distributed 
among different locations |50W53| . 

In this paper, we restrict our investigation to the fol¬ 
lowing cases: We study position-momentum measure¬ 
ments or their optical equivalent, quadrature phase am¬ 
plitude measurements. First, we analyze the steering of a 
two-mode squeezed state, to understand the decoherence 
effects of coupling to a thermal reservoir. We examine 
two properties of the decoherence: the effect of simple 
losses (damping); and the effect of thermal noise. Sec¬ 
ond, we consider a special case of EPR-steering, that 
relates to a macroscopic/mesoscopic superposition state 
- a “Schrodinger cat” state. In fact, the EPR-steering 
paradox can be used to signify the Schrodinger cat su¬ 
perposition [54] . We examine the effect of decoherence 
as caused by interaction with the environment, on this 
signature. 

An interesting new feature evident for EPR-steering 
is that the nonlocality can manifest asymmetrically with 
respect to the observers [551150] . In our results, we focus 
on the asymmetry of the decoherence effects on the EPR- 
steering. This is valuable to understand, not only from 
the fundamental perspective of testing quantum mechan¬ 
ics, but from the point of view of the potential applica¬ 
tions of EPR-steering, which include cryptography EEK 
166] and no-cloning teleportation [65]. One-sided device¬ 
independent cryptography has been proposed using EPR 
steering inequalities [53]. In recent papers, it was shown 
how the asymmetric effects of loss on EPR-steering could 
affect the optimal location of teleportation stations m, 
or which entanglement criterion should be used to faith¬ 
fully verify entanglement in cases of untrusted devices 

[SMB]- 


II. DECOHERENCE OF STEERING FOR THE 
TWO-MODE SQUEEZED STATE 

Consider two quantum harmonic oscillators, with bo¬ 
son operators a and b, and depicted in the Figure 1 as 
systems A and B. Such oscillators can be coupled, so 
that EPR-steering correlations are induced between the 
two oscillators. 


One of the nicest examples of such EPR steering has 
been realised experimentally in optics m as the two¬ 
mode squeezed state |69| . Here, each mode of the light 
field is modeled as a quantum harmonic oscillator. An ex¬ 
ample of a coupling between the oscillators that generates 
this type of state is called parametric down conversion, 
and can be described by the interaction Hamiltonian in 
an appropriate rotating frame m 

H = inE(ab — a 1 b^) (1) 

Similar two-mode squeezed states can be generated by 
impinging two single-mode squeezed states into the two 
different input ports of an optical beam splitter im 

The two-mode squeezed state system is fundamentally 
interesting, as it can be quite accurately realised in op¬ 
tics, but also it gives a reasonable approximation to the 
effects we expect to see in many other physical realisa¬ 
tions of steering, such as in opto-mechanics Enumi, 
atomic ensembles |19H21| . and BEC |22fJ30] , The study of 
the effect of the thermal environment on the EPR steer¬ 
ing for a two-mode squeezed state will therefore give us 
insight into a broad set of physical scenarios. The two¬ 
mode quantum correlation effects were noticed originally 
in the context of photon number correlations between the 
two beams, which gives rise to noise reduction mm- 

A. The two-mode squeezing Hamiltonian 

We begin by reviewing the solution for the EPR cor¬ 
relations given by the Hamiltonian of Eq. 0 . The 
EPR solutions were originally derived in [16] . We de¬ 
fine the quadrature phase amplitudes Xa, Pa, Xb and 
Pb for each mode: Xa = a + a\PA = (a — cd)/« and 
Xb = b + , Pb = {b — b^)/i. This choice of scaling for 

the amplitudes will imply the Heisenberg uncertainty re¬ 
lations AXaXPa > 1 and AXbAPb > 1. We note that 
depending on the physical system modelled by the Hamil¬ 
tonian, the amplitudes can also correspond to scaled po¬ 
sition and momentum observables. We now ask what 
correlations are generated between the quadrature ampli¬ 
tudes, after an interaction time r between the two modes, 
at sites denoted by A and B. 

On examining ([I]), the resulting coupled equations for 
a and b can be readily solved, to give 

a( T ) = r/a(0) - ^/(r/ 2 - l)6 f (0) 

b(r) = r,6(0) (2) 

where ry = coshr. We define the squeezing parameter as 
r = The quadrature phase amplitudes are given by 

Xa (t) = coshrA^ (0) — sinhrXg (0) 

Pa (t) = coshrPA (0) + sinh rPs (0) 

Xb{j) = coshr As (0) — sinhrA^ (0) 

Pb ( t ) = cosh rP B (0) + sinh tPa (0) (3) 

This enables us to calculate correlations such as (A^), 
( XaXb ) after a time t, assuming initial uncorrelated 
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vacuum states for all modes. We can then evaluate 
the variance of X A — g x X B , which will be denoted 

A 2 (X A — g x X B ) and which equals (^(Xa — g x X B ) 2 ^ 
when we assume vacuum inputs at r = 0. We have intro¬ 
duced constants g x , which can be chosen to minimise the 
variance relative to the Heisenberg quantum noise level 
of the amplitudes X'a and Pa- We find this value of g x 
by standard procedures mm-- 

±'A‘ 2 (X A -g x X B )=0, ^ 9x = ^^l ( 4 ) 

For the system described by the Hamiltonian of Eq. (1), 
we find the optimal value is g x = — tanh2r. The optimal 
variance is given by: 

A 2 {X A - g x X B ) = (X\) {X f*f (5) 

which for the parametric system ( 1 ) becomes 
A 2 (Xa - g x X B ) = cos j i 2r . Similarly we evaluate 
the variance of Pa + g P PB , which we will denote as 
A 2 (P A + 3 p Pb): The variance minimizes when 

—A- (P A + g P P B ) = 0, g P = - ( 

This optimal value of g p for system (1) is g p = — tanh2r. 
The optimal variance is given by 

A 2 (P A + g P P B ) = < P \) - {P *p^ ( 6 ) 

which becomes A 2 (P A + g p P B ) = cos ^ 2r for Eq. (1). 

B. EPR steering correlations 

The criterion for the EPR paradox introduced in m is 
also a criterion for EPR steering mg. This EPR steering 
criterion is defined as the square root of the variance 
product: 

EPRa\b(9x, g P ) = A (X A - g x X B ) A (P A + g P PB) (7) 
We observe EPR steering (of system A by B) whenever 

EPR A \B(g x ,9 P ) < 1 ( 8 ) 

The ideal correlations created by the parametric down 
conversion Hamiltonian interaction (1) are obtained in 
the limit of r —> oo and correspond to the EPR variances 
becoming zero, and hence EPR —> 0. 

The EPR steering condition © by its very definition 
will negate all local hidden state (LHS) models of the 
type El 3 

(X e A Xt) = (^4) = / P(X)(X e A ) x (X^ B ) x dX (9) 


These LHS models are similar to the local hidden vari¬ 
able models introduced by Bell ia mg. Here A represents 
hidden variable parameters and P( A) is the probability 
distribution for these parameters. The P( A) is indepen¬ 
dent of the choice of measurement (6 and </>), which is 
made after the generation and separation of the subsys¬ 
tems A and B. In this model, the {X A ) x is the average 
value for the result X A , given the hidden variable state 
specified by A. The (X B ) x is defined similarly. For the 
LHS model however, there is the additional asymmetrical 
constraint that the local hidden variable moments (such 
as (X A ) x ) for system A be consistent with measurements 
of some local observables (for example position and mo¬ 
mentum) at A, and is thus able to be described as arising 
from a local quantum density operator p A . 

We can define the minimum value of EPR after opti¬ 
mising g x and g p as: EPR A \b = min{EPR A \ B (g x ,g p )}. 
Thus, for the Hamiltonian of Eq. 0 : 

EPR «° = < 10 > 

as derived originally in m- Ideal EPR correlations are 
obtained as r —> oo, in which case EPR A \b ~t 0 , and we 
say the “steering is perfect”. 


C. Reservoir coupling 

We now examine the effect on EPR steering if the sys¬ 
tems A and B are independently coupled to heat bath 
reservoirs (Figure 1). We assume that the parametric 
interaction H given by (1) that generates a two-mode 
squeezed state is turned off at the time t = 0 , and the 
system left to decay. The two-mode squeezed state is a 
so-called Gaussian state, meaning that its characteristic 
function is Gaussian m- Within the constraint of two¬ 
mode Gaussian states and measurements, the condition 
(I 8 | will provide a necessary and sufficient test of steering 
[3j ! 4j. This makes the criterion valuable, for understand¬ 
ing the effects of decoherence in Gaussian systems. 

Thus, we consider a system prepared in a two-mode 
squeezed state at time t = 0. In principle, the modes a 
and b can be spatially separated. We consider the cou¬ 
pling of each mode a and b to independent heat baths 
(reservoirs) with thermal occupation numbers n a and 
respectively. The solutions after coupling to the reser¬ 
voir are straightforward to evaluate using the operator 
Langevin equations that describe the evolution of the 
mode operators CUES]: 

a = -7aa+V / 27^fA(t) 

b = -^bh+ x /2j b rb(t) ( 11 ) 

Here q Q and 7 b describe the decay rates (losses) that 
are induced by the reservoirs. The quantum reser¬ 
voir operators T(t) have nonzero correlations given by 

<rt (t)T a (t')} = n a S(t - t') and <r a (i)rt (if)) = (n a + 
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1)5 (t — t') where the numbers n a and rib give the level of 
thermal occupation of the reservoirs. Solutions are 


a(t) = e-^aM+V^U \ a (H) dt' ( 12 ) 


From this, we can calculate the moments at a later time, 
in terms of the initial moments: 


(a!(t)a(t)) = e- 2 ^ t (a t ( 0 )o( 0 )) + n a { 1 - e" 2 ^ 4 ) 
(a(t)b(t)) = e~ ( ^ + 7 b)t (a(0)b( 0)) (13) 


The solution for (b^(t)b(t)) is obtained from that for 
(a' (t)a(t)), but exchanging the letters a with b. The ini¬ 
tial moments are given by the solutions found in Section 


II.A. The final moments after reservoir coupling are 

(a\t)a(t)) = n a (l — e~ 2t7a ) + e _2t7a sinh 2 (r) 

( a{t)b(t )) = _ e ~ha.+ib)t coshr sinhr (14) 

The covariance matrix V is defined as: 


V = 


(X\) (X A P A ) (X A X B ) (X A P B ) 
(PaX a ) (Pi) (P A X B ) {P A P B ) 
(X B X A ) (.X B P A ) (A 2 ) (. X B P B ) 

(P B X A ) (P B P A ) (P B X B ) (Pi) 


(15) 


Hence we find V 14 = V 23 = 0 , V r2 = V 34 = i , Vu = V 22 , 
V 33 = V 44 , V 24 = — V 13 where 


V 11 = e~ 2lat cosh 2 r + (l — e _27 “ 4 ) (1 + 2 n a ) 

V 33 = e~ 2lbt cosh 2 r + (l — e~ 2lbt ) (1 + 2 rib) 

V 24 = e~^ ,a+ ' yb ^ t sinh 2 r (16) 


Applying the results of the Section II. A, we calculate that 


EPR A \b — 


cosh2r[e 2 ' Vat (l — e 2jbt ) (1 + 2 nb) + e 2lbi (l — e 27at ) (1 + 2n a )] + e + (1 4- 2n a ) (1 + 2nj,) (l — e 2 ^ bt ) (l — e 


e 2 ^ bt cosh 2r + (1 — e 2 T'b t ) (1 + 2 rib) 


(17) 




Figure 2. The decoherence of two-mode EPR steering with 
no thermal excitation ( n a = rib = 0). A sudden-death type of 
decoherence is noted when it is the system “doing the steering” 
that is is coupled to the reservoir. The steering signatures 
EPR A \b and EPR b \a are plotted versus t' = 'ybt, for various 
values of squeeze parameter r. EPR steering is obtained for 
the Gaussian system iff EPR < 1 and the strongest steering is 
when EPR -A 0. Left graph: Reservoir coupling to system 
B only ( 7 a = 0). Solid lines correspond to EPR A \ B while 
dashed lines correspond to EPR B \ A . Curves are (bottom to 
top on the far left) r = 2 (black), r = 1 (red), r = 0.5 (blue). 
Right graph: Plots of EPR A \ B for symmetric coupling to 
the reservoirs: = 7 b- In this case, EPR A | S = EPR B \ A . 


D. Decoherence of EPR steering with no thermal 
noise 


of loss to be asymmetrical with respect to the two sys¬ 
tems A and B. 


1. Damping the steering system: Steering sudden-death 


Let us suppose that the reservoirs are coupled asym¬ 
metrically, so that 7 a —t 0 i.e. only the mode b is coupled 
to a reservoir. In this case, we find for the steering of A 
(mode a) 


EP R A \ b — 


e 2lht + cosh 2 r [l — e 27fc4 ] 
1 + e~ 2lbt (cosh 2 r — 1 ) 


(18) 


The plots of Figure 2 (left graph) indicate that where 
the losses are entirely on the steering mode B, the de¬ 
coherence of EPR steering is substantial. After a time 
given by = m, steering (as measured by EPR A \ B ) 
in that direction is lost. The loss is inevitable, regardless 
of how much steering is present in the initial two-mode 
quantum system, as shown by the results for the differ¬ 
ent values of squeeze parameter r: We also note that the 
cut-off time for steering is independent of the amount of 
steering r in the initial two-mode system. The behavior 
observed here is analogous to the “entanglement sudden- 
death”, that has been observed for the decoherence of 
entanglement [761178] , 


We first assume negligible thermal noises: we put 
n a = rib = 0 in the solution given by Eq. (17). The 
decoherence effect on the system by the reservoirs mani¬ 
fests as losses, parametrised by y a and 7 b- We will study 
the effect on the steering parameters, showing the effect 


2. Damping the steered system 

However, the sudden-death effect is not apparent when 
the “steered” system is lossy. We again assume y Q —> 0 


— 2 7a* 
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and evaluate the steering of the damped system, by the 
undamped system: 


EPRb\a — 


e 2lht + cosh 2 r [l — e 27i,t ] 
cosh 2 r 


(19) 


The plots of Figure 2 (left graph) reveal the nature of the 
decoherence when only the steered system is damped. 
Here, we see that there is much less sensitivity of the 
steering to the losses. In fact, while the EPR steering is 
certainly diminished by the dissipation, there is no cut¬ 
off, or sudden death, but rather steering is still possible 
for arbitrary times, t -A oo. We remark that the in¬ 
creased sensitivity to losses affecting the steering system, 
as compared to the steered system, has been noted in 
earlier papers CUES, SB, 591 . 


3. Symmetric decoherence 


We next consider a reservoir with symmetric damping 
7 a = lb- 


e " 47bt + (1 - e~ 2lbt ) 2 
A ^ B 1 + e -27i,t (cosh 2 r — 1 ) 

2 cosh 2 r \e~ 2lbt — e -476 *] 
1 + e -276 * (cosh 2 r — 1 ) 


( 20 ) 


The results of Figure 2 (right graph) show, as might be 
expected, the effect is dominated by the fact that the 
steering system is damped. Here steering (as measured 
by EPRa\b) is also lost at 7 b t = 



Figure 3. The decoherence of two-mode EPR. steering with 
a thermal reservoir coupled to system B only. Here y a = 0, 
r = 1. We plot the steering signature versus t! = y b t. The 
thermal noise creates a sudden-death type decoherence effect 
when it is the system that is “being steered” that is coupled to 
the reservoir. Solid lines correspond to EPR A \ B while dashed 
lines correspond to EPR B \ A . EPR-steering is obtained when 
EPR < 1 and the strongest steering is for EPR —► 0. We 
plot different values of n b : from bottom to top (for each line 
type), n b = 1 (blue), n b = 5 (red), n b = 10 (black). 



E. Decoherence of EPR steering with thermal noise 

Now we consider the more complex interaction where 
there is additional thermal noise for each reservoir. Here 
we use the full expressions given by Eq. (17). Our re¬ 
sults in Figures 3 and 4 reveal several features. We note 
that the loss of steering is rapid and complete (“sudden- 
death”) if the thermal noise is placed on the system that 
is being steered. This effect has been noticed in previous 
studies of thermal steering and is especially important for 
asymmetrical systems such as found in opto-mechanics 

(2D EH 13111351. ' The effect is more significant as the ther¬ 
mal noise increases, but is reduced for higher r (Figure 

4). 

F. Comparing with the decoherence of 
entanglement 

The effects of the reservoir on the steering are asym¬ 
metrical. This will not be true for entanglement: En¬ 
tanglement is symmetrically defined, with respect to the 
two different systems involved, and the numerical value 
we assign to the entanglement will be unchanged under 


Figure 4. The decoherence of two-mode EPR steering with 
thermal reservoirs. EPR.-steering is obtained when EPR < 1 
and the strongest steering is for EPR -A 0. Left graph: 
Reservoirs coupled to both systems A and B symmetrically. 
We plot EPR A \ b vs t! = y b t. Here y a = 7 b, n a = n b , r = 
1. We plot different values of n b , given from bottom to top 
(for each line type) , n b = 1 (blue), n b = 5 (red), n b = 10 
(black). Solid lines correspond to r = 1 while dashed-dotted 
lines correspond to r = 2. Right graph: Thermal reservoir 
coupling to system A and cold reservoir coupling to B (y a = 
76 , n a = 1,5,10). Solid lines correspond to EPR A \b while 
dashed lines correspond to EPR b \ A . 


the exchange A <—> B. The relationships between en¬ 
tanglement and EPR steering for Gaussian states have 
been recently studied EUSQ]. For Gaussian states, it is 
possible to give a quantification of steering ElfflESHSg, 
which is an asymmetrical quantum correlation, similar to 
quantum discord |811(55] . 

We may evaluate the entanglement via the parameter 
Ent = A(X a - gX B )A(P A + gP B )/{ 1 + g 2 } EH IB]. 
The inequality Ent < 1 is a necessary and sufficient con¬ 
dition for entanglement of two-mode Gaussian systems 
for the case we examine here where A (X A — gXs) = 
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Figure 5. The decoherence of the entanglement of the two- 
mode squeezed state under the action of a reservoir. The 
entanglement signature for Gaussian systems Ent is plotted 
versus scaled time t ', under different sorts of reservoir cou¬ 
pling. Entanglement is obtained iff Ent < 1 and stronger 
entanglement corresponds to Ent —» 0. Left graph: Here 
we take all thermal noise zero ( n a — nt = 0). Solid curves 
show the symmetric case = 76 - Dashed-dotted curves are 
for the asymmetric case, for a reservoir coupling to system B 
only (y a = 0 ) when t' = •ybt', or, for a reservoir coupling to 
system A only ( 7 ^ = 0) when t! = 'y a t. We plot different val¬ 
ues of r, given from bottom to top (for each line type): r = 2 
(black), r = 1 (red), r = 0.5 (blue). Right graph: Here we 
study thermal reservoirs. Solid curves show the symmetric 
case n a = nb and 7 a = 7 b- Dashed-dotted curves are for the 
asymmetric case, for a reservoir coupling to system B only 
( 7 a = 0) when t 1 = 'ybt. We plot different values of thermal 
noise, given from bottom to top (for each line type) n 6 = 1 
(blue), nb = 5 (red), nb = 10 (black). 

A (Pa + 9Pb)- The condition for optimally chosen g has 
been shown equivalent to the Peres-Simon condition in 
that case m■ Thus, in this section, we consider the 
entanglement parameter 

Mx A - g x B) Hr A + 9 P B ) 

1 + 9 2 

First we find the value of g that minimizes the param¬ 
eter Ent. We evaluate d Qg t , in this case we know that 
A (X A - gX B ) = A (P A + gP B ) hence: 

dEnt d A 2 (X A — gX B ) 

dg dg 1 + g 2 

_ ~2 (X A X B ) + 2 g ((A 2 ) - (X 2 )) + 2 g 2 (X A X B ) 

(1 + 5 2 ) 2 

( 22 ) 

Next = 0, 

dg ’ 

O - (. X A X B ) + g ((X 2 ) - (X 2 )) + g 2 (X A X B ) = 0 
Hence the value of g that minimizes Ent is ES3: 

- «X 2 ) - (X 2 )) + V /((X 2 )-(X 2 )) 2 +4(X A X B ) 2 
9 ~ 2 (X A X B ) 

( 23 ) 


We note that if (X^) = (X B ), then the value of g that 
minimises the expression is g = 1, but where we have 
asymmetric reservoir effects, like different reservoir cou¬ 
plings or thermal noises, the optimal value will be differ¬ 
ent to 1. Here we will use from Section II.A that: 

(X 2 ) = e _27at cosh2r + (l — e~ 2lat ) (1 + 2 n a ) 
(X A X B ) = -e" ( 7 “ + 76 )t sinh2r 

(X|) = e~ 2lbt cosh2r + (l - e- 27tt ) (1 + 2 n b ) 

(24) 

Figure 5 plots the entanglement Ent for various reser¬ 
voir couplings. We see from Figure 5 (left graph) that 
where there is no thermal noise, the entanglement de¬ 
cays steadily, but is never destroyed completely. This 
effect was noted in [55] [57] . The larger values of r corre¬ 
spond to larger amounts of entanglement for the initial 
state before decoherence, and we note that while the de¬ 
cay is sharper for higher r, a larger initial amount of EPR 
entanglement will ensure larger EPR entanglement for all 
later times. The Figure 5 (right graph) shows that when 
the reservoirs are thermally excited, the entanglement is 
totally destroyed (in sudden-death fashion) after a finite 
time. This time is shortened, by increasing the amount 
of thermal excitation of the reservoir (on either system). 

III. DECOHERENCE OF THE STEERABILITY 
OF A S-CAT STATE 

Traditionally, most studies of decoherence have cen¬ 
tred around the Schrodinger cat |45H15] . We noticed 
in the study of the decoherence of steering for the two¬ 
mode squeezed state that for stronger EPR effects, the 
decay was brought about more sharply, to give overall 
decoherence times of a similar order (Figure 2). This is 
consistent with the overall intuition about decoherence 
in quantum mechanics, that it acts to destroy the more 
extreme (larger) effects more quickly so that they are 
not observed in nature. This has been studied for the 
Schrodinger cat example, where one considers a system 
initially prepared in superposition of two states macro- 
scopically distinguishable e.g. in phase space | l.'S'll!')] . 
The coupling to a reservoir induces a decay of the super¬ 
position. Where the separation in phase space increases, 
calculation and experiments show that the decay rate 
will increase. Thus, there is an explanation of the transi¬ 
tion from microscopic to macroscopic quantum mechan¬ 
ics. Here, motivated by this, we examine the decoherence 
of the EPR steering of a Schrodinger cat state. 


A. Steering signature for a Schrodinger cat 

Consider the state: 

W = ^{\-a) A \t) B +P s \a) A \i) B } (25) 
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-4-3-2-10 1 2 3 4 


Pa 


Figure 6. The conditional distributions associated with an 
EPR steering signature for the Schrodinger cat. Left Graph: 
Conditional probability for X A given that the result of a mea¬ 
surement a z is +1 (red solid curve is for a = 2, green dash- 
dotted curve is for a = 3) or —1 (blue long dashed curve is 
for a = 2, black short dashed curve is for a = 3). Right 
Graph: Conditional probability for Pa given that the result 
of the measurement a x is +1 (red solid curve) or —1 (blue 
long dashed curve). Here a = 2. 


which is the entangled Schrodinger cat state. We select 
9 = 7 t/2 and a to be real. Here, \a) is the coherent state 
for mode a of system A and | f), | J,) are the eigenstates 
of the Pauli spin a z of system B. This state is similar to 
that described in Schrodinger’s original gedanken exper¬ 
iment j39| . where a microscopic system becomes entan¬ 
gled with a macroscopic one, the “cat”. The spin-mode 
entangled state given by Eq. (251 has been the subject 
of several experiments |43l . We first explain how one can 
signify the Schrodinger cat, using EPR steering. 

We will derive the probability distributions P(X A ) (or 
P(Pa)) f° r the system A, given that a measurement of 
the Pauli spin a z (or a x ) at B has been carried out to 
give a result +1 or — 1. We define the scaled amplitudes 
for position and momentum by a = (x + ip) where a 
is the boson operator, and note that for real position 
and momentum, c = V2. In terms of the quadrature 
phase amplitudes X A /b and Pa/b defined in Section II, 
however, we chose the scaling c = 2. We will distinguish 
the two cases by using lower and upper case, respectively, 
and note we have dropped the use of the operator hats 
where the meaning is clear, or to denote the outcomes of 
the measurements. As expected from direct examination 
of the state (25), the distributions P(x|+1) and P(x| —1) 
for the result x at A given the outcomes +1 or —1 for a z 
at B, are the two Gaussian hills (Figure 6 , left graph): 
Specifically, 


/ 

\ / 2 'i 

21 

„ x 2 „ 9 xa 

( X 

II 

t-H 

-H 

II 

N 

b 

-exp 

- 2 — - 2a 2 T 4— 

\ 

/ \7T y 

c 

c 


(26) 


Next, we derive the conditional distributions for P(P A ) 
given that a measurement of Pauli spin a x at B yields 
the result +1 or —1. To evaluate this, we rewrite the 
state in terms of the basis states | and | |) x for the 



B 


Figure 7. The steering of a Schrodinger cat. We have pro¬ 
posed a signature for the cat, based on EPR steering param¬ 
eter EPRa\b- Here, measurements made on the microscopic 
spin system B “steer” the macroscopic system (the cat) A. 
The “cat” in this case is modelled as the superposition in¬ 
volving two coherent states, |a) and — a). The value of a 
determines the size of the Schrodinger cat. We consider that 
the “cat” being macroscopic is coupled to a thermal reservoir. 
We examine the decay of the steering signature, with respect 
to the parameters of the reservoir coupling. 


spin a x . 

W = \i\ ~a) A + i\a) A )\ t )x,b 

+ ^{i\a) A -\-a) A )\l)x,B (27) 

At this point, we note the for other choices of 9 , normal¬ 
isation factors are more complicated. We find that the 
probability for the momentum given the results for spin 
a x is 


(p a x = ±l) 


2\ 2 1 o p2 f ( P a 

-e~ 2 P* (^1 ± sin ^47 


(28) 


This allows us to calculate the conditional variances 
that give us a signature for EPR steering. An EPR steer¬ 
ing of the “cat” is observed when 


EPR = Var{X A a z )Var(P A o x ) < 1 


(29) 


where here 

Var{X A a z ) = P z {+l)Var{X A +l)+P z (-l)Var{X A -1) 


denotes the conditional variance for X A averaged over the 
outcomes of a z . We have used the notation Var(X A ±1) 


to mean the variance of P(X A ±1), which is the proba¬ 


bility distribution for X A given the outcome ±1 for the 
measurement a z , respectively. The Var(P A \& x ) is de¬ 
fined similarly, but for outcomes of a x . The right-side 
bound is the quantum noise limit, as determined by the 
Heisenberg uncertainty relation which in this case for the 
























choice c = 2 is AXbAPb > 1. We note that the deriva¬ 
tion of the inequality as a signature of EPR steering has 
been given in the Refs. ®\mm- 

The EPR quantity ( [29] ) can be evaluated: We see that 
for the two Gaussian hills, the variance is at the quan¬ 


tum noise level (Figure 6 , left graph): Var{X a 


cr z ) = 1 . 


On the other hand, the variance associated with the mo¬ 
mentum distributions is reduced below 1 (Figure 6 , right 
graph), implying that EPRa\b < 1- We can calculate 
the variance from the conditional distributions. Alterna¬ 
tively, noting that the collapsed state for system B given 
a measurement a x is the superposition |(| — a) A + i\ot) a) 
(for result + 1 ) or a — I — &)a ) (for result — 1 ), it is 

easy to use the methods and results of the next section, 
to find that 


Var(PA cr x ) = 1 — 4a 2 e 


(30) 


The steering inequality (291 has been suggested in 


Ref. [54], as a way to realise an EPR paradox with the 
Schrodinger cat state. The steering cannot be obtained 
if the system is in the mixture of states, | — ot) a\ t)s and 
| a) a | that allows a classical dead or alive descrip¬ 

tion. This case is especially interesting, because it focuses 
on the steering of a mesoscopic system (the Schrodinger 
cat) [H]. If we assume Local Realism is valid, then it is 
the local state of the mesoscopic system that is shown 
to be inconsistent with the quantum mechanics (refer 
to the LHS model of the Section II.B). This contrasts 
with EPR’s original argument, which showed the incon¬ 
sistency for a microscopic system [7]. 


B. Decoherence of the Schrodinger cat with a heat 

bath 


It will prove useful to next study the interaction of the 
single mode “Schrodinger cat” state [15H15] 

\^ = ±{\-a)+z\a)} (31) 

with a reservoir. This was analysed by many authors, 
including Yurke and Stoler |4 5114 01 . We consider that 
the single mode, prepared in the Schrodinger cat state, 
is then coupled to a thermal heat bath with dissipation. 
Using the solutions from Section II.C, we can calculate 
the moments at a later time t, in terms of the initial 
moments: 

(a^(t)a(i)) = e~ 2lt (a^(O)a(O)) + n (l — e ~ 2lat ) 
(as(t)a(f)) = e~ 2lt (a(O)a(O)) 

Denoting the variances using shorthand notation, by 
A 2 P = (A P) 2 and A 2 A' = (AX) 2 , we find 

A 2 P = 1 + 2n (l - e" 27t ) - 4a 2 e~ 27 t e ~ 4a2 (32) 



0 0.05 


0.1 

t' 


0.15 0.2 


Figure 8 . The decoherence of the Schrodinger cat EPR- 
steering signature EPR for the system depicted in Fig. 7. 
Here, a thermal reservoir is coupled to the spin system B. We 
plot the steering signature versus t! = 7 1. Here we consider 
different sizes of the “cat” (solid line lower to top, a = 0.5, 
1 and 1.5). For each a, we plot different values of n: n = 0 
(solid), n = 1 (dashed-dotted). 


We see that the variance for P can be reduced below the 
quantum limit (given by 1). The reduction of the vari¬ 
ance of P below the quantum limit is in itself a signature 
for the “cat” state. We see this as follows. If we take 
a large, and assume no thermal noise and 7 = 0 , then 
the probability distribution P(x) associated with each of 
the “dead” and “alive” states is a Gaussian hill, with vari¬ 
ance 1. If the system were to actually be in any kind of 
mixture of these quantum states, then the variance for P 
could not drop below 1 because of the uncertainty rela¬ 
tion (and the fact that mixing states cannot decrease the 
variance) [54l [ 88 ]. The observation of a reduced variance 
for P can thus occur for a superposition but not for a 
mixture of the two quantum states that possess a distri¬ 
bution P(x) given by the Gaussian hills. The details are 
not studied further here however, since our objective is 
to study the decoherence of steering. 


C. Decoherence of the steering of the Schrodinger 

cat 


Motivated by this, we now examine the decoherence of 
the EPR steering of a Schrodinger cat state, as given by 
the signature Eq. (291. We consider the two-mode case 


where A is a harmonic oscillator coupled to a heat bath 
from time t = 0 and the second system is the spin B (no 
heat bath coupling), as depicted in Figure 7. The system 
is prepared in the entangled Schrodinger cat state 


\ip) =-^{\-a} A \f)B+i\a)A\l)B} (33) 


We can rewrite this state in terms of the basis states | f)x 
and | 4_) x for the spin a x as shown in Section III.A. It is 
straightforward to calculate the moments (aa z ), ( a 2 a z ), 
etc. at the time t = 0. The moments at a later time t 
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(after the interaction with the heat bath has been turned 
on) can be evaluated using the solution (12 1 of Section 


II.C. We note that since the reservoir interaction does not 
involve the spins, we have a z (t) = er 2 ( 0 ), <r x (t) = <r x ( 0 ). 
Therefore, we can calculate the moments at the time t, 
in terms of the moments at the initial time t = 0 . 

We next consider the moments of the distributions 
at the time t for X A and Pa, conditional on getting 
the result either +1 or —1 for the Pauli spin measure¬ 
ment on B. These moments can be readily calculated. 
We outline the method. We use the notation, for ex¬ 
ample, (Xa 1 ) z to denote the moment (Xa) of system 
A , conditioned on the result +1 for spin <r 2 at B. Let 
us suppose then that we obtain the outcomes Xa at 
A, and a z at B. We define the measurable probabil¬ 
ity P(Xa,(J z ) for the joint outcomes. We see that: 


P(X A 

P(X A 


<t 2 ) = p (j’) • It follows that P(X A ,±1) = 
±1)P*(±1) = \P(X A ±1). Hence, (X A a z ) = 


P(X A ,1)X A - P(X A ,-1)X A = P(X A 1 )P z (1)A a - 

P(X A \-1)P Z (-1)X A = H^l 1 ) - <*a|- 1». Follow¬ 
ing this procedure, we see that: 

(X A a z ) = - 2 {(X A \l) z - (X A \-l) z ) 

(Xa) = \{(Xa\i) z + (X a \- 1) z ) (34) 

and similarly for the moments involving X\. This 
allows us to solve for the conditional moments using 
the relations such as (X A l) z = (X A cr z ) + (X A ) and 

(X A —1 ) z = —(X A cr z ) + (X A )- We also define the 
same relations, but replacing X A with P A and a z with 
a x . The final solutions for the conditional moments are: 
(X A \±l) z = Tlae-*, (A 2 |± 1} 2 = 1 + 2n (l - e~ 2 ^) + 

4a 2 e~ 2yt , (P A ±1)* = P 2 ae^e " 2 " 2 and (P\ ±l) x = 
1 + 2n (l — e _27< ). Hence we solve for the average vari¬ 
ance Var(X A <J Z ) of A conditioned on the measurement 
of spin a z at B: We obtain 


Var(X A c z ) = 1 + 2n(l — e 27t ) 


(35) 


Similarly, we can solve for variance Var(P A <J X ) of P A 
conditioned on spin a x at B, to obtain 


Var(P A &x) = l + 2 n (l — e 27t ) — 4a 2 e 27t e 


—4a 


(36) 


To show how the EPR steering signature decoheres with 
time, we evaluate 

EPR = Var(X A \a z )Var(P A \a x ) (37) 

Here, we use the notation defined in Section III.A. 


We find that with n = 0 (no thermal noise), the value 
of EPR is identical to that of Var(P A cr x ), which is also 

identical to the variance A 2 P given by Eq. (32). This 
value is plotted in Figure 8 (solid lines). By the above 
argument, the signature of the S-cat is the drop below 1 
of EPR. We note that as a increases, the EPR value 
tends to 1 , the signature therefore becoming more sen¬ 
sitive to decoherence. This reduction in variance of P 
as a increases is directly associated with the interference 
fringes in the distribution P(P A ) |45fl48| . These fringes 
become finer as a increases, as evident from the function 
given by Eq. (28). With thermal noise, we see that the 
value of EPR increases: a much greater sensitivity to 
decoherence is apparent, as illustrated in Figure 8 by the 
dashed-dotted curves. 


IV. CONCLUSION 


The directional properties of EPR steering are signifi¬ 
cant when it comes to understanding the decoherence of 
EPR-steering. We have shown how when two systems 
are coupled independently to a reservoir, the steering is 
asymmetrically affected. If we consider the steering of 
a system A by measurements made on the system B, 
then the steering is sensitive to the reservoir coupling to 
system B. This is intuitively not surprising, given the 
Local Hidden State (LHS) definition of EPR steering as 
a nonlocality. In terms of the LHS description that is 
to be negated in order to confirm such steering by B 
m, it is the hidden states of system B that are like 
those considered by Bell. The sensitivity of the steering 
to the losses on the system B has been studied experi¬ 
mentally, in the context of the EPR paradox [B 6 : , %7|. 
There, it was known that the EPR criterion could not be 
achieved, with 50% or more losses on the steering sys¬ 
tem A. Recent work considers this result in terms of 
the monogamy properties of EPR steering [55]. The sen¬ 
sitivity of the steering signatures to losses (which may 
represent an eavesdropper on the channel) has poten¬ 
tially important implications for quantum cryptography 

mm- 

The behaviour with respect to the thermal noise ap¬ 
pears almost reversed. This has implications for de¬ 
tecting steering of Schrodinger cats which are coupled 
to hot reservoirs. For the Schrodinger cat example, we 
explained how the steering acts as a signature of the 
Schrodinger cat. Adding thermal noise to the system B 
significantly affects the amount of EPR steering of sys¬ 
tem A. There have been recent theoretical studies of the 
steering of a mechanical oscillator, which we liken to a 
“cat”, by measurements made on an optical pulse fMll55] . 
In those studies, the thermal noise of the oscillator was 
shown to have a quite dramatic effect on the amount 
of steering possible. By comparison, the entanglement 
between the thermal oscillator and the pulse was quite 
robust. Similar results have been obtained for Bose Ein- 






























stein condensates m- These results are consistent with Dalton for many helpful suggestions, 
the simple descriptions of decoherence given in this pa¬ 
per. 
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